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ABSTRACT 

The  necking  process  of  an  axisymmetric  tension  specimen  made 
of  elastic-plastic,  strain-hardening  material  is  analyzed  by  a  generalized 
J 2* flow  theory  for  large  deformations.  The  governing  equations  are 
solved  in  a  Kantorovich  type  approach  based  on  a  variational  principle. 

The  effects  of  both  geometric  nonlinearities  resulting  from  large  de¬ 
formations  and  the  physical  nonlinearities  arising  from  plastic  material 
behavior  are  considered.  Numerical  results  have  been  found  for  the 
stress  and  deformation  histories  in  the  specimen  up  to  a  50%  reduction 
of  the  neck  radius.  The  shape  of  the  neck,  for  the  first  time,  is  rationally 
calculated. 
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I.  Introduction 


The  conventional  tension  test  in  iiHnally  conHidcred  to  provide  the 
basic  leiiHile  properties  of  materials.  The  lest  does  provide  knowledge 
of  the  intrinsic  stress-strain  relation  up  to  a  point.  However,  after 
nec  king  occurs  in  a  tension  specimen,  the  distributions  of  stress  and 
strain  become  very  complicated  and  are  not  known.  The  information 
obtainable  from  the  test  thereafter  is  only  some  kind  of  average  stress- 
average  strain  relation.  A  typical  tension  specimen  after  necking  occurs 
is  shown  schematically  in  Fig.  1. 

In  order  to  provide  a  better  understanding  of  the  true  stress- strain 
relation,  Bridgman  |1J  studied  the  stress  distribution  at  the  narrowest 
section  (section  A-B  in  Fig.  I )  of  a  necked-down  axisymmetric  tension 
specimen  and  obtained  a  semi-empirical  solution.  His  analysis  was 
based  on  the  assumption  that  the  axial  strain  across  the  narrowest  section 
was  uniform.  Davidenkov  and  Spiridonova  |<i],  by  an  analysis  similar  to 
Bridgman's,  obtained  another  result  for  the  stress  distribution  at  the 
neck.  In  both  of  tnese  analyses  the  displacements  and  hence  the  shapes 
of  the  necked  profiles  of  the  tension  specimen  were  not  calculated.  In 
fact,  the  determination  of  the  stress  distribution  had  to  rely  on  measured 
values  of  the  neck  diameter  d  and  radius  of  curvature  p  (see  Fig.  1). 

These  solutions,  therefore,  were  not  entirely  rational.  Thomason  [3] 
has  calculated  the  profiles  of  the  axisymmelrically  necking  tension 
specimen,  based  on  a  "maximum  unloading  rate"  hypothesis.  But  no 
solution  for  the  stress  distribution  was  included  in  his  analysis.  Recently, 
Segal  [4]  has  studied  the  plastic  flow  in  the  neighborhood  of  the  neck  of  an 
ideally  plastic  tension  specimen.  He  made  some  calculations  both  for  the 
stress  and  the  displacement  fields.  However,  the  analysis  did  not  provide 
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FIG  1  NECKED-  DOWN  TENSION 
SPECIMEN 
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a  complete  solution  for  the  stress  and  the  strain  distributions.  These 
earlier  investigations  evidently  do  not  provide  an  adequate  analytical 
solution  to  the  problem  of  necking  in  strain-hardening  bars. 

In  a  rigorous  analysis,  it  is  necessary  to  find  the  displacement 
and  the  stress  in  the  whole  specimen,  not  just  at  the  neck;  and  for  the 
complete  loading  history,  not  just  at  a  particular  deformation  stage. 

The  effects  of  plastic  unloading  must  be  considered;  also,  the  large 
changes  in  geometry  that  occur  during  necking  must  be  taken  into  account. 

The  purpose  of  the  present  analysis  is  to  calculate  rationally  the 
deformation  (and  hence  the  necking  shape)  and  the  stress  distribution  in 
a  tension  specimen  throughout  its  loading  history.  The  tension  speci¬ 
men  considered  is  of  finite  length  and  initially  has  some  imperfection 
on  the  lateral  surfaces  in  such  a  way  that  the  -adius  at  the  middle  is 
slightly  smaller  than  the  radius  at  the  ends.  The  specimen  is  subjected 
to  a  uniform  axial  displacement  at  the  ends  and  is  elongated  in  such  a 
way  that  the  ends  remain  shear  free.  The  specimen  material  is  assumed 
to  be  isotropic,  elastic -plastic,  and  stra'n-hardening.  A  variational 
equation  of  equilibrium  for  large  deformation,  in  conjunction  with  the 
approximation  method  due  to  Kantorovich  [5],  is  used  to  solve  for  the 
field  variables.  A  generalized  ^"flow  theory  for  elastic -plastic 
materials  with  large  deformation  effects  is  used  in  the  analysis.  Thus, 
both  the  effects  of  the  geometric  nonlinearities  resulting  from  large 
deformations  and  the  physical  nonlinearities  arising  from  plastic  material 
behavior  are  included  in  this  analysis. 

II.  Field  Equations 

A  set  of  field  equations  for  analyzing  arbitrary  deformations  of 
elastic -plastic  materials  will  be  presented.  The  equations  involved  are 
written  in  a  general  form  with  the  use  of  tensor  notation. 
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2.  1.  Equilibrium 

The  notation  and  equation*  of  equilibrium  used  in  the  present 
analysis  follow  those  of  [6],  The  pertinent  results  of  [6]  will  be  briefly 
reviewed. 

Let  a  material  point  of  a  solid  body  be  identified  by  the  coordi- 

12  3  — 

nates  (f  ,  f  £  )  and  let  the  vectors  £  and  £  denote  respecitvely,  its 

position  when  the  body  is  undeformed  and  deformed.  The  covariant 

*JL 

base  vectors  of  the  undeformed  and  deformed  body  are  =  — r  ,  and 

1  »€ 

?  .  =  — r.  respectively,  and  the  covariant  components  of  the  metric 

1  as1 

tensors  are 


g. .  =  e  . 
*1J 

G. .  =  ¥ . 

ij  ~i 

*  ij 

with  inverses  g1^  and  G1^.  The  contravariant  base  vectors  are  <£  -  g1^. 

J 

and  =  G^^j. 

One  defines  the  contravariant  tensor  components  u1  of  displacement 
u^  =  £  -  £  by  the  relation 

u  =  u*  e  .  (1 ) 

The  Lagrangian  strain  tensor  is  defined  by 

T)..  =  i(G..  -  g..)  =  i(V.u.  +  V.u.)  +  ig  ♦(*•  uS)(V.ut)  (2) 

*ij  2'  ij  6ij'  2'  j  i  i  j'  26st'  i  ''  j  '  ' 

where  V  denotes  covariant  differentiation  in  the  undeformed  body.  The 

stress  tensor  jr  in  the  deformed  body  may  be  written  symbolically  in 

terms  of  its  contravariant  components  <r1^  as 


<r  =  (r1^  e  .  e  . 


(3) 
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Consider  an  arbitrary  volume  V  of  the  undeformed  body  bounded 
by  a  surface  S;  after  deformation  these  become  V  and  S,  respectively. 
Let  n  be  the  unit  normal  to  S  and  define  n.  by  the  relation 

~  i  ' 


H  =  niJL  * 

If  body  forces  are  neglected,  the  equation  of  force  equilibrium  governing 
the  stress  rate  and  displacement  rate  given  in  [6]  is 

V.«lir  +  <iij?.ur  +  qij?.ur)  =  0  (4) 

where  q^  =  —  <r^  is  the  Kirchhoff  stress;  here  G  and  g  are  the  deter¬ 

minants  of  G^j  and  g„,  respectively,  and  where  (*)  denote  differentiation 
with  respect  to  a  monotomcally  increasing  quantity  such  as  time  t. 

•  *  j 

If  surface  traction  rates  ^  =  T  per  unit  original  area  are 
prescribed,  then 


(qir  +  +  q^n.  =  Tr  (5) 

on  the  boundary.  Prescribed  displacement  rates  v,  =  vl£.  are  satis¬ 
fied  by  letting  ul  =  v^. 

A  variational  equation  governing  u1  and  q1^,  which  can  be  deduced 
from  the  classical  principle  of  virtual  work,  is 


J  {«l1^C»?ij  +  |q1^[gst(viuS)(^.ut)]}dV  =  y  T^u.dS  (6) 

V  S 


for  all  states  q1^,  ii1,  T1  satisfying  equilibrium,  and  all  displacement 
variations  6u.,  with 

i’ 


«hij  =  +  (^jut)(V.6u!')]} 


Either  the  equilibrium  equation  (4)  or  the  variational  equation  of  equili¬ 
brium  (6)  can  be  used,  in  principle,  for  analyzing  plasticity  problems 
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involving  large  deformations.  In  the  present  analysis,  the  variational 
i-quation  (6)  will  be  used,  in  conjunction  with  a  Kantorovich  approxima¬ 
tion  method,  to  solve  for  the  field  variables. 

he  field  equations  are  still  incomplete  and  must  be  augmented 
by  constitutive  relations  between  stress  and  strain--or,  more  precisely, 
between  stress  rates  and  strain  rates. 

2,2.  Constitutive  Equation 

There  have  been  several  forms  of  the  constitutive  equations  for 
elastic -plastic  materials  involving  large  deformations,  such  ns  those 
suggested  by  Green  and  Naghdi  [7],  Lee  and  Liu  [8],  Cameron  and 
Scorgie  [9],  Lee  [10],  and  Willis  [11].  However,  in  the  present  analysis 
of  the  large  deformation  in  a  tension  specimen,  the  new  constitutive 
equation  given  below  seems  to  be  more  convenient  to  use. 

The  materials  considered  are  isotropic,  elastic -plastic  and  strain¬ 
hardening.  For  simplicity,  neither  a  Bauschinger  effect  nor  the  influence 

of  temperature  is  included.  It  is  assumed  that  the  actual  strain  rate  rj. . 

.  E 

can  be  written  as  the  sum  of  the  elastic  strain  rate  rj. .  and  the  plastic 

p 

strain  rate  17.^.  The  elastic  stress -strain  relation  as  suggested  by 
Budiansky  [12]  is 


i  ik 

where  <r.  =  a  G^,  is  the  mixed  tensor  component  of  the  stress  tensor, 

E  is  Young's  modulus  and  v  denotes  Poisson's  ratio.  This  is  a  genera¬ 
lized  Hooke's  law  which  includes  large  deformation  effects.  Equation  (7) 
is  not  integrable  and  is  not  a  rigorous  constitutive  relation  derived  from 
the  theory  of  elasticity  for  arbitrary  deformation.  But  it  will  not  result 
in  significant  errors  for  the  present  analysis,  because  for  infinitesimal 
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strain  Equation  (7)  is  the  exact  statement  of  Hooke's  law;  and  for  large 

deformation  in  the  present  necking  problem,  the  elastic  strain  only 

plays  a  minor  role  and  is  dominated  by  the  plastic  strain. 

The  plastic  stress-strain  relation  is  a  generalization  of  the 

flow  theory  to  include  large  deformation  effects.  Let  us  begin  with 

the  flow  theory  for  infinitesimal  strains.  Referred  to  the  Cartesian 

coordinates,  the  stress  and  strain  tensor  components  are  <r. .  and  e... 

*  iJ  iJ’ 

the  stress  deviator  is  s..  =  cr. .  -  —6 .  .<r.  .  ,  and  the  effective  stress  is 

ij  ij  3  ij  kk’ 

3  — 

defined  by  <r  =  f4s..s..l2.  The  yield  function  is  F  =  <r  -  c,  where  c 
e  1 2  ij  ijJ  7  e  * 

is  the  initial  yield  stress  or  the  maximum  value  of  tr^  which  has  occur¬ 
red  in  the  stress  history,  and  the  flow  rule  states  that  the  plastic  strain 
rate  is 

ij  2'Et  E  <rfi  ij  e  e  dt'e' 


if  <r  <  c, 
e 


or  tr  =  c  and  ir  ^  0 
e  e 

where  E(  is  the  tangent  modulus,  which  is  the  slope  of  the  uniaxial 
stress-strain  curve. 

Next,  let  large  deformation  effects  be  included.  Referred  to 
general  coordinates,  the  contravariant  components  of  stress  are  <r^, 
and  the  contravariant  components  of  the  stress  deviator  tensor  £  are 

-  icij<rklG 

s  0-  3u  0-  ukl  . 

The  mixed  tensor  components  Sj  are  defined  as 

i  i  1  *i  kl— 

J  J  3  j  kl 


i 
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win*  ro  6!  G,kGk  .s  llu*  Knmoikor  do  It  a  and  eklGkj  is  the  mean  normal 
citrons.  Tlio  effective  stress  is 


it  (4M-- 

o  <2  j  i 


l4(«r,j<rklG..G.1  -  *U‘j<rklG..G  )H  . 
£  tk  jl  3  U  kl 


jl  J  ij 

and  tho  \iold  function  is  still  F  -  e  -  o.  The  cciicrali/od  flow  rulo 

o 


(«) 


st.itos  that 


%  -  E«-'k,0ki0ij  •  7Vkl°ki'^  “*o  -  *•  '"d J  0 


<9)a 


i)**  -  0 
»J 


if  *  <  o  . 
o 


or  9  -  c  and  fr  *♦  0 

e  f 


<9)b 


where.  as  shown  later,  the  tangent  modulus  E(  is  now  defined  as  the 
slope  of  the  true  »tress*natural  strain  curve  in  uniaxial  tension.  Note 
that  the  plastic  strain  gi\en  by  (<*)a  satisfies  the  plastic  incompressi* 
bilitv  condition.  This  can  be  seen  by  examining  the  rate  of  change  of 
volume  which  is 

■  &<•  ?  dV>  ■  r.  ?  °lJVv 


-  f  c‘J(i,o  *  i’!jwv  ■ 

The  change  of  volume  contributed  by  the  plastic  strain  rate  is*  2  G^h^tiV, 
which  by  (°)a.  vanishes. 

Combining  equations  (?)  and  (9)  gives  the  full  constitutive  equation 
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'ij  *  E1!*'  *  v,(Ctk*i  *  cjk*!'>  ‘  %*£  I  “  \  “■  *n<1  *„  •*  0 


.  3 .  I  I  w  kl_  _  I  -  kl_  % 

mz)i9  ckicij  -  iV  C^% 


<10 1* 


^  *  Eli"  ♦ ♦  vN  ■  *Vk>  ‘  • 

or  ff  =  «  *nd  fr  0 

A  <1 

no»b 

Physically,  Equation  (I0)a  applies  in  the  plastic  loading  region  and  flfjfb 
applies  in  the  plastic  unloading  or  elastic  region.  The  strain  rate  given 
by  (I0)a  depends  not  only  on  the  ci  rrent  stress  rates  hut  also  on  the 
whole  stress  and  deformation  histories. 

In  order  to  define  E(  precisely,  let  us  now  consider  uniaxial  tension. 
The  stress-strain  relation  then  reduces  to 


*>n  •iicu»!i*ir  -eI»"giigh  ■? 


if  v  •  c  and  k  +  *t 

0  fS 

i  Ilia 


I  -  .  1 

*11  E°u*i 


if  <  C  , 


or  e#  ■  c  and  k*  «  0  flllb 

where  the  effective  stress  is  f*  *Cj  j  •  rj,  by  Equation  (6).  Since  r*  * 
is  the  tensor  component  referred  to  the  base  vector  £|,  which  has  the 
magnitude  Cj  j  .  the  corresponding  physical  component  of  stress  is 


9  «  e1 1  ‘  C, 

XX  1  I 


cll  "l  •  \ 


Therefore,  Equations  (1 1  )a  and  (I  l)b  become 


ft 


k 

xx 


M2, a 


The  natural  strain  is  defined  by 


r)  =  ln(j--) 


ln(  1 


(11) 


where  u%  is  the  displacement  along  the  coordinate  x;  !q  and  I  are  the 
original  and  the  current  length  of  a  cylindrical  bar  subjected  to  uniaxial 
tension.  The  rate  of  change  of  the  natural  strain  is 


1 


04) 


By  comparing  (12)  and  (14)f  we  obtain  the  following  simple  relations  for 
the  uniaxial  tension  case: 


* *  “  •««  • « *Bd  *«« ' 0 

"■^KK  “  \x<C* 

or  9  ■  c  and  ¥  <0 

xx  xx 

Now  the  tangent  modulus  C(  is  clearly  seen  to  be  the  slope  of  the  true 
stress-natural  strain  curve  as  shown  in  Fig.  2,  and  E  is  the  linear  re¬ 
lationship  between  q  and  when  plastic  unloading  occurs  (as  well  as 
during  loading  in  the  absence  of  plastic  flow).  The  constitutive  equation 
(10)  is  used  under  the  condition  that  the  *xx  -  q  relation  (or  equivalently 
E  and  E()  in  uniaxial  tension  is  known. 

The  constitutive  equation  is  complete,  but  its  inversion  (the  stress 
rate  in  terms  of  the  strain  rates)  is  still  to  be  found.  This  inversion 
must  be  used  in  conjunction  with  the  variational  equation  (6),  because 


i 
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we  wish  to  wrtte  *lj  in  terms  of  rj^,  end  successively  write  in 
terms  of  ii1  end  solve  for  ii*.  The  inversion  of  (10)  is  derived  in  the 


Appendix;  the  result  is 


#lj(  »  .1.  >1. 


"  £ >•  Out 

..li  E  _ik_.il.  3  Et  E  kl 


*  II  ♦  ,{ff  -  Iv!  cUo>tI|»kl  * 


if  e#  »  c  end  s^hjj  **  0  (15)* 


‘1J  •  irfrt  °ikcJ\i  ♦  (i  ♦  »!ff  -  i.»G‘Jck\i 

•  <.ikcJ't>kl  ♦  •,kC,ll|tt|l  U  ».  <  c  . 

or  ■  c  end  s^Sj^  **  0  (I5)b 

The  field  equations  ere  now  complete  for  analysing  lerge  deforme* 
tion  problems.  Note  thet  the  Kirchhoff  stress  (|^  eppeers  in  the  equa* 
tions  of  equilibrium  wherees  the  constitutive  equetions  involve  In 
the  present  analysis,  however,  qIJ  cen  be  epproximeted  by  For 
smell  deformetions,  C  *  g.  so  thet  q^  ■  e^  —  r'K  In  the  ne:king 
problem  lerge  deformetions  occur,  but  since  the  plestic  strein.  which 
is  incompressible,  dominetes  the  elestic  strein,  the  volume  chenge 

C  dV 

ceused  by  the  totel  strein  is  very  smell.  Then,  since  —  •  m  ly 
it  is  still  true  thet  q1^  ■  *lK  Therefore,  q^  will  be  repleced  by 
in  the  equetions  of  equilibrium  for  simplicity. 

Finally,  by  the  use  of  the  constitutive  equetion  (15)  we  cen  obtein 


another  version  of  the  variational  statement  thet  turns  out  to  be  convenient 
for  numerical  analysis.  It  cen  be  easily  verified  thet 
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*%  * 

both  for  plastic  loading  and  unloading.  Therefore,  if  1  is  defined  by 

1  *  i  £  «*\  ♦  -  £  T'u,  dS  ,16, 

v  5t 

where  Sj  ia  the  part  of  the  original  surface  on  which  traction  rates  are 
prescribed,  then 
6l  «  0 

is  the  same  as  the  variational  equation  (6). 

111.  A  Perfect  Bar  Under  Uniaxial  Tension 

Before  we  go  into  the  necking  problem,  let  us  study  the  uniform 
deformation  of  a  perfect  circularly  cylindrical  bar  under  uniaxial  tension. 
The  bar  has  an  initial  length  2L  and  an  initial  radius  Rq  aa  shown  in 
Fig.  3.  The  ends  of  the  specimen  are  subjected  to  a  prescribed  uniform 
axial  displacement  U^.  The  radial  displacement  at  the  ends  is  unre¬ 
strained  so  that  the  ends  remain  shear  free  during  the  elongation.  The 
following  is  the  solution  for  the  uniform  deformation  of  the  perfect  bar; 
this  solution  will  be  used  for  further  comparison. 

The  physical  axial  displacement  of  the  specimen  is 
ut  *  a  •  a 

u, 

where  a  *  is  the  engineering  axial  strain.  The  natural  axial  strain 
is  n  *  ln(l  ♦  s).  For  a  given  engineering  strain  a,  the  uniaxial  stress 
in  the  specimen  can  be  calculated  by  the  uniaxial  stress-strain  relation. 

In  this  paper,  we  will  use  the  Ram  berg -Osgood  law  (13)  to  relate  the 
uniaxial  true  stress  »  and  the  natural  strain  r)  as  follows: 


(17) 
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where  n  is  the  strain-hardening  exponent,  and  K  is  a  material  constant. 

3  o-  j  _n 

As  shown  in  [13],  K  is  equal  to  -y(-|r)  ,  where  <r  is  a  secant  yield 

strength  at  which  the  secant  modulus  E  =  (~)  =  .  7E.  In  the  present 

s  ™ * 

n 

analysis,  it  will  be  convenient  to  use  the  dimensionless  stress  quan- 
<r 

z  z 

tity  (— — ).  Rearranging  equation  (17)  gives  the  following  relation  between 

<r  s 
z  z 

( — — )  and  a  given  value  of  a : 
s 

=  .  (17)a 

s  s 

The  radial  displacement  is  u^  =  br,  where  the  engineering  radial  strain 
b  can  be  found  by  the  use  of  the  constitutive  equation  (10)a.  Let  be 
the  component  of  the  Lagrangian  strain  tensor  in  the  r  direction  and  let 
<r33  be  the  only  non- vanishing  contravariant  component  of  the  stress 
tensor  (i.e.  consistant  with  the  cylindrical  coordinate  system  £*  =  r, 

^  =  0,  and  £3  =  z);  then  rj ^  ^  =  b  +  -|b2  and  =  <r3  =  =  <r  . 

One  of  the  constitutive  relations  of  Eq.  (10)a  reduces  to 


<r  fr 


<r  <r  ,  & 
Bw  zz.n-1,  zz 


*11  ■  '  Ti  n0ii'^)'Ti>n'  If5* 


where  the  metric  tensor  =  (1  +  b)  .  Integration  of  it  gives 

b  =  e  s  8  -  1  (18) 

The  radius  of  the  specimen  after  deformation  is  R  =  Rq(  1  +  b).  Next, 

the  total  load  in  dimensionless  form  is 

<r 


T  =  1-^X1  +b)2 

s 


JTR?<t 
0  s 


,  zz  . 

=  (  — )e 

s 


(19) 
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<T 


Not» 


o  tliat  tlu>  values  of  b,  (  Z—  )  and  ( - 4 — )  depend  on  three  material 

(r8  - 


JlRnir 

0  K 


parameters  only;  namely,  the  strain-hardening  exponent  n,  the  ratio 


<r 


(-=r)  and  Poisson's  ratio  v. 

<r 

For  the  particular  values  n  8,  (-|r)  0.  0074  (eonsistant  with 

the  values  E  -  <25.8  x  10^  psi.  and  <r  -  .  186  x  10^  psi.  for  Cr-Ni  steel), 

s 

1  T  —  T 

and  v  -j ,  curves  relating  a  and j —  *  and  ^  and 1 —  are  plotted  in 

itR^o-  7TR  n0" 

0  s  0  s 

Fig.  4  and  Fig.  5.  A  very  simple  approximation  for  the  maximum  value 

of  the  total  load  (found  by  neglecting  elastic  strains)  is 

1 

_  „  c r  ■— 

1  -2  'max  1 7  *  E  ^ 

"Vs 


and  the  corresponding  values  for  the  stress  and  engineering  strains  are 

1 


.  22 .  f3  .  s n 

(  o-  T=Tmax  “  ^  7  n  E  ^ 


(a)T=Tmax  =  *  “  1 


_  _1_ 

(b),_  „  =  e  2n  -  1 

'  ;T=Tmax 


<r 

g 

This  approximation  should  hold  fairly  well  for  n  >  2  and  n(-=r)  «  1. 

For  n  =  8,  for  example,  ^a)j-Xmax  *s  *  and  the  corresponding 
natural  axial  strain  is  (n)j_jmax  =  •  125. 

IV.  Formulation  of  the  Necking  Problem 

The  tension  specimen  considered  for  the  present  necking  problem 
is  an  imperfect  bar  as  shown  in  Fig.  6,  of  which  the  lateral  surfaces 
deviate  slightly  from  the  straight  surfaces  of  the  perfect  bar  of  Fig.  3. 
The  initial  length  of  this  imperfect  bar  is  again  2L,  but  its  initial  radius 


ENGINEERING  AXIAL  STRAIN  FOR  PERFECT 


PERFECT  BAR  ( R0*1) 
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ft 


% 


2 

'7'  '  < 

la  R(i)  R^(i  *  <c  |,  where  <  (assumed  small)  Is  the  maximum 
deviation  between  the  surfaces  of  the  perfect  and  imperfect  bar  at 
i'0,  and  a  la  a  poaitive  decay  parameter. 

The  apecimen  la  aaaumed  to  be  axisymmetric  about  the  x*axi« 
and  to  be  aymmetric  with  respect  to  the  middle  plane  a  *  0.  The  oniy 
non* vanishing  physical  components  of  displacement  are  uf  and  u#, 
perpendicular  and  parallel  to  the  s-axis,  respectively.  The  non- 
vanishing  phyairal  components  of  atreaa  are  rfr,  e^.  m  and  ffi. 

The  atreaa  and  displacement  components  are  functions  of  r  and  a  and 
the  time-like  variable  t;  they  are  independent  of  •  by  axial  symmetry. 
The  boundary  conditions  imposed  at  the  ends  are  the  same  as  those 
imposed  on  the  perfect  bars  namely,  at  the  ends,  u^fr,  tL,t)  ■  tu^(t). 
Also,  the  radial  displacement  at  the  ends  is  unrestrained  so  that  the 
ends  remain  shear  free  during  the  elongation. 

Referred  to  the  cylindrical  coordinates  (r,  t,  a),  the  metric  tensor 
g(j  has  components 

*11  *12  *11  1 

*|j‘  *21  *22  *21  *  0 

*11  *12  *11_  0 

The  contravariant  and  covariant  components  of  the  displacement  ere 


U  •  U|  •  uf 


u  •  u  j  •  ut 


and  the  stress  tensor  has  contravariant  components 

.11  *  _I1 


9 

« 


0  e 
_22 


11 


.11 


0 
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Thc  rates  of  the  Legranittan  strain  tensor  components  are 

du  eu  du  du 

'll  '  "  *7r> 7?  *  7T  7T* 


hu  «  (r  ♦  ur»or 

#u  tfu  iu  du 

*»•"  t1  *r 


uoi 


,  do  do  ,  do  du  ,  du  di»  du  do 

hu  'ii  ‘i"  *  TT'ir  *  7r'»r  ‘j'lrir*  »r»r’ 

The  covariant  components  of  the  metric  tensor  of  the  deformed  bod)  are 

*ur  2  *“/  2 

Cn  •«'  ♦  TT'  *'?f' 


CU  *  <r  *  V 


*u.  2  ®“r  2 

°J5-"  *»^> *«tF» 


(21) 
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with  the  Im 
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«o.i°ii-0.i' 


i 

C  *  3T” 


’22 


122) 


G33. 
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(ciicn  *  °i  i] 


G1  3  «  C31  i 


-g 
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The  non-vanishing  contravariant  components  of  the  stress  deviator, 

II  22  J1  .  13 

s  .  s  ,  s  and  s 


are  given  by 


J 


%,J  ,‘j  -  ic'-,<#,,C|I  ♦  *UCU  ♦  *J,Gn  ♦  2*  *  *0, 


Tin  •>§!••*  Iiw«»  «trv»k  t* 


-  I  *n  »Mr-*  22.2_2  »Mr-*  *  •  »  2#*2 

•,.  )C||*|%  I  cw  ♦  N  I  C}i  ♦  <(k  I  C|  j 

. .  1 1.2_  ~  .  ,  II  >>,.2  .  .  II  I  ,. 

*  il'  1  ci icji  ^  %  °hM'  %  cnci> 

♦  4*,,%l,ClJCJJ|)-  U4 

For  th<  |>rv»vnt  problem,  l«*t  >i«  «  onttnuc  to  •<%•?  IK*  RkiiJn  ri'*Okfcv/l 
liiw  (17)  <•%  thv  unikklal  true  klrvkk*n«lural  strain  D.UtlMi.  IK*;  taiiv  M 
motlif'ii<k  l«y  definition  l» 


I  n_  1  1  J  ,Vn-l 

F  ‘  V  '  E  ♦  E  ln{Ti 


t  v  % 

In  thv  following  wv  will  uav  thv  convenient  dlmenalonk-aa  atrva*  quantity 

JJ  %IJ  \ 

(— )  and  *1*0  the  dimenatonlca*  quantitive  I  — )  and  I  — |.  Suoetit  mon 
•  a  % 

of  (25)  into  (15)  give* 

•j;'  •Trbi<f  Hc"‘c>,iV|-p(^M^M^.kk,i 


*  (1  ♦  v)(l  -  2v)  ^  *C  C  ^hl  **  e#  ,C  fy*.!  1  e#  ,C  °>,1^ 


where 


27n(^)n“J 

■—  ■  *  If  9  «  c  and  »^o  .  *  0  f27)a 

%  n-l  *  »J 

26(1  ♦  v)  ♦  I2n(-*)n 


P(~) 


If  e^  <  c.  or  «  c  and 

ij.  -  r 

•  %  1 


I27ib 
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The  constitutive  rsUtlons  in  matrix  form  are 
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C.2 

CI3 

CH 
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S 

^22 

9 

cn 

C22 

C23 

C24 

"22 

S 

frJ3 

9 

S 

1 

cj. 

C32 

C33 

C34 

*33 

fr,J 

e 

C4. 

C42 

C4  3 

C44 

2h,3 

where  the  coefficients  are 


II  ,  9 


it 

9  .^11 


e„  ■  rh«f  <o">2  -  .i-»2p«^»l  - 


cM  ■  rh<f*ir^t  *  <^>2*V»i  -  *<£«22 
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c«  *  JiTTTi'f'i0' *  rhz  <C‘J>2  *  i«V^»2P<^H 

a  s  s 
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s  s  e 
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II  22  9 
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» 


■24* 


C24  *  C42  * 
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^"g1*  •  (7-)(V-»p<r 

S 
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11  i* 

7-)Gij 

♦(^)G11) 

s 

s 

_  ,,  ..  22  33  9 

s 

s  s  s 

ifiirf 

,,  1*  22  13  0 

s 
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s 

It  tt  »  .33  » 
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|(—^)G33  ♦  (^)G13) 
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With  the  use  of  the  constitutive  equation  (21),  the  equation  (Is)  of 
the  variational  principle  becomes 

L  *(•> 

^  *  21  J0  J0  tCl  l^fl  f  f  *  4C44i,13  f  2ClZi,llhiZ 

*  2CI  3**l  I  **33  *  40 14^1 1  **13  *  *  4C24I,22,|1  3 


II  du  ,  du  ,  22  , 

♦  «CMn„b„  ♦  ♦  (-5T>  I  ♦  <77»4,> 


.11  *1.  ,  til.  ,  .11  til.  til.  til.  tii. 

♦ «77'H»r>  - <ir»  i * wir ♦ sr st *>rdrd*  • 

(30) 


The  surface  integral  of  (16)  vanishes,  since  the  prescribed  tractions 
vanish  on  all  boundaries.  Equation  (30),  with  the  use  of  the  strain* 
displacement  relation  (20)  can  be  written  in  the  form 


25 


*  if  \  \  [QiMdrittc  terms  in  u  ,  u  ,  and  ihstr 

»  J0  J  0  r  * 

derivative  sjrdrds 

Tne  variational  principle  fll  *  0,  with  the  I  of  (51),  will  enable  one  to 


solve  for  t>r  and  &  . 


Finally,  the  geometric  boundary  conditions  that  must  be  imposed 
on  the  displacements  are 

i»g(r,L,t)  -  UL(t)  (3 

at  the  end,  and 


ut(r,0,t)  ■  0 

at  the  middle,  by  symmetry. 


The  governing  equations  of  the  present  problem  were  approximately 
solved  by  the  Kantorovich  approach  (5)  based  on  the  variational  principle 
6l  ■  0,  where  1  is  given  by  Equation  (31).  In  this  method  the  displace* 
ments  were  assumed  in  the  form 


ur  ■  rf(t,t)  ♦  r  n(s,t)  (34) 

u#  «*<*»*>♦  r2k(*,t)  (35) 

where  f,  h,  g,  and  k  are  functions  of  the  space  variable  s  as  well  as  the 
time  t.  The  axisymmetric  conditions  require  that  equation  (34)  contain 
only  odd  powers  of  r  and  that  equation  (35)  contain  only  even  powers  of  r. 
To  analyse  the  necking  process,  we  must  calculate  the  stress 


and  the  displacement  u^  as  functions  of  r,  s  and  t.  The  details  of  this 
calculation  will  be  shown  later;  the  following  paragraph  outlines  briefly 
the  computation  steps  involved  in  the  procedure. 


« 


-£b~ 


i 


Assuming  knowledge  of  ihe  histories  of  <r  J  and  it  at  time  t, 
calculate  the  displacement  rate  (u.)t  by  the  variational  equation  for  a 
prescribed  axial  displacement  rate  at  the  end.  Then  calculate 

and  by  the  strain*displacement  relations  and  the  constitutive 

equations,  respectively.  Next,  extrapolate  to  find  (<r13),,  »'  d 

<u.). .  ...  and  then  use  this  knowledge  to  calculate  (ii.).  .  ...  Repeat 

I  itAK  1  (tAl 

the  same  procedure  for  successive  small  time  intervals  and  thus  obtain 

the  stress  and  displacement  histories  for  a  finite  time  span. 

By  this  procedure,  the  displacement  functions  f,  g,  h,  and  k  of 

Equations  (34)  and  (35)  are  assumed  to  be  known  at  time  t,  and  it  is  their 
•  •  • 

rates  f,  g,  h.  and  k  that  are  the  unknowns  to  be  found.  Substitution  of 
the  assumed  displacements  (34)  and  (35)  into  the  equations  (20)  for  the 
rates  of  the  Lagrangian  strain  tensor  components  gives 

i|||  ■  (I  +  f  +  3r2h)(f  +  3r2h)  ♦  4r2kk 


h22  =  (1  +  f  ♦  r2h)(f  ♦  r2h) 

*»33  8  T\  +  r2hz,(iz  +  r\>  +  (l  *  8z  *  r\)(*z  +  T\' 

r»i  3  s  h3i  s  2  r({z  +  +  3r2^  +  \  r<*  ♦  *  ♦  3r2h)(£x  +  r2hz) 

♦  rk(gz  +  r2kz)  +  r(l  ♦  gz  ♦  r2kz)it  (36) 

With  the  use  of  these  relations,  the  expression  (31)  for  the  integral  I 
of  the  variational  principle  becomes 


rt FJ5**»** 
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j-* 

T-  f„  lDll^+DlA+D13*  +  D14»'z+D15ie*+D16ik  +  D17iki 
s  *y0 


•  ^  ••  ••  ••  •  •  •• 

+  D„r  +  D„f  h  +  D,,f  h  +  D,  A  g  +D,,f  ktD„f  k 
22  z  23  z  24  z  z  25  zBz  26  z  27  z  z 


+  D33h2  +  D^hh,  +  D^hg,  +  D,Ahk  +  D„hk 


-Will  T  .  CI1K  I  _  /  i 

34  z  35  Bz  36 


37  z 


+  D44hz  +  D45hz®z  +  D46hzk  +  D47hzkz 


*2 


+  D55gz  +  D56gzk  +  D57gzk2 


+  D66k  +  D67kkz 


+  D77k^]dz 


(37) 


where  the  coefficients  DQp  are  definite  integrals  with  respect  to  r.  For 
example 

RU) 

Dn  =  2?r  ^  {Cn(l  +  f  +  3r2h)2  +  C22(l  +  f  +  r2h)2 

+  2C12(1  +  f  +  3  r2h)(  1  +  f  +  r2h)  +  2C14r(l  +  f  +  3r2h)(fz  +  r2hz) 

+  2C24r(l  +  f  +  r2h)(fz  +  r2hz)  +  C44r2(fz  +  r2hz)2 

11  ^  22  2,  . 

+  <r  +  cr  r  /  rdr 

R(z) 

D12  =  271  J  {2(C13  +  C44>r2(l  +  f  +  3r2h)(fz  +  r2hz) 

+  2C23r2(l  +  f  +  r2h)(fz  +  r2hz)  +  2Cur(l  +  f  +  3r2h)2 


+  2C24r(l  +  f  +  r2h)(l  +  f  +  3r2h)  +  2C^r3(f,  +  r2hj2 


34*  '  z 


1  ^ 

+  2r<r  }  rdr 
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The  rest  of  these  28  coefficients  are  calculated  in  a  similar  fashion, 

but  the  results  are  not  shown  here.  All  the  D  a  are,  of  course,  func- 

ap  7 

tions  of  z. 

The  condition  6l  =  0  can  then  be  imposed  in  order  to  solve  for  f, 

•  • 

g,  h,  and  k,  as  functions  of  z.  Let  the  integral  domain  of  (37);  namely, 

the  plane  bounded  by  z  =  0,  L  and  r  =  0,  R(z),  be  divided  into  m  x  n 

meshes,  with  m  +  1  mesh  points  in  the  r  direction  and  n  +  1  mesh 

points  in  the  z  direction  (see  Fig.  7).  Numerical  calculation  at  time  t 

starts  by  the  computation  of  G..,  G1^,  S1^,  <r  and  C..  by  means  of 

®  ^  J 

Equations  (21),  (22),  (23),  (24),  and  (29)  at  each  mesh  point.  In  com¬ 
puting  C„,  we  have  to  know  whether  a  mesh  point  is  governed  by  the 

plastic  loading  conditions  or  by  the  plastic  unloading  conditions  in  order 
<r 

to  calculate  P(  — )  by  (27)a  or  (27)b,  respectively.  However,  before  the 
s 

solutions  for  u^  are  found,  these  plastic  loading  or  unloading  conditions 

cannot  be  checked,  and  hence  the  position  of  the  boundary  between  the 

plastic  loading  region  and  the  plastic  unloading  region  is  not  known.  In 

the  present  analysis,  this  boundary  is  determined  by  an  iterative  process 

which  will  be  described  later.  For  now  we  just  assume  a  trial  position 

tr 

of  the  boundary  in  the  specimen  and  calculate  P(~)  either  by  (27)a  or 

s 

(27)b.  After  C„  is  evaluated  at  each  mesh  point,  we  integrate  numeri¬ 
cally  along  the  r  direction  by  Simpson's  rule  to  get  the  DQp  of  (37).  Then 

th 

let  the  discrete  values  of  f,  g,  h,  k  and  at  z^  (i  mesh  point  along 
the  z-axis)  be  f.,  g^,  h^,  k.  and  (DQp)„  By  the  use  of  the  trapezoidal 
rule,  Equation  (37)  can  be  written  approximately  in  the  summation  form 


n 

I  =  ^  L[Quadratic  te 


•  • 


rms  inf.,  Bm,  h.,  h.+1>  k.,  k.+1, 

(38) 


<Vi  and  <Vi+ll 
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The  variational  principle  6l  =  0,  with  the  I  of  (38),  implies  that 

JL-.o 

3f. 

i 


JL  =  o 

agi 


91 

8h. 


=  0 


91 

8k. 


=  0 


(39) 


for  i  =  2  to  n.  The  natural  boundary  conditions  of  the  variational 
principle  are 
81 


=  0 


8f 


1 


81 

8h, 


=  0 


(40) 


at  the  middle  z  =  0,  and 


81 


=  0 


8f. 


n+1 


81 


=  0 


8h 


(41) 


n+1 


at  the  end  z  =  L.  In  addition,  the  geometric  boundary  conditions  (32)  and 
(33)  become 

*1  -o 

kj  =  0 


fen+1  =  UL 


kn+l  =  0 


(42) 


T 
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Equations  (39),  (40),  (41)  and  (42)  constitute  a  system  of  4(n  ♦  1) 

•  •  • 

linear  equations  which  contain  4(n  4  1)  unknowns  f.#  ft.,  and  ^  (for 

i  =  1  to  n  +  1).  These  equations  were  solved  by  Potter's  method  (see 

•  • 

Ref.  14)  with  the  aid  of  an  IBM  360/65  computer.  After  f.t  g.,  )r  and 

k.  are  found,  the  displacement  rates  u  and  u  can  be  calculated  by 
i  r  z 

means  of  (34)  and  (35). 

The  calculated  displacement  rates  u^  and  u^  arc  not  necessarily 
the  correct  solution,  because  they  were  obtained  on  the  basis  of  an 
assumed  position  for  the  boundary  between  the  plastic  loading  and  plastic 
unloading  region.  So  after  we  calculate  rj.j  from  it,  we  must  check  the 
plastic  loading  conditions  — namely,  o-^  =  c  (c  is  the  maximum  value  of  «r 
which  has  occurred  in  the  stress  history)  and  rj.jS1^  **  0— at  each  mesh 
point  in  the  assumed  plastic  loading  region.  If  these  conditions  are  not 
satisfied,  we  redefine  the  mesh  point  assumed  to  be  plastically  loading 
to  plastically  unloading.  Similarly,  we  check  the  plastic  unloading  condi¬ 
tions —  namely,  <re  <  c,  or  =  c  and  r).jS1J  **  0— for  each  mesh  point  in 
the  assumed  plastic  unloading  reg'on  and  make  the  required  adjustment 
if  necessary.  Thus  we  obtain  a  new  position  of  the  boundary.  Then  we 
calculate  u.  on  the  basis  of  the  new  boundary  and  do  the  checking  again. 

If  this  checking  is  successful,  the  last  displacement  rate  u.  found  is 

taken  as  the  right  solution.  From  these  values  for  u.,  the  strain  rate 

1 

rj. .  and  the  stress  rate  (  — )  are  calculated  by  (36)  and  (28),  successively. 

s 

The  stress  and  the  displacement  at  the  time  i  t  At  are  estimated 
by  a  simple  parabolic  extrapolation  of  the  form 

Vit  =  Ut  +  lI4t  '  K-/J4* 


assuming  a  knowledge  of  the  various  rates  at  the  time  t  -  at  (initially 
at  t  *  0,  it  suffices  to  use  *  uQAt  if  At  is  sufficiently  small)  Fo r 
convenience,  the  ttme*like  variable  t  was  identified  as  the  prescribed 
end  displacement  UJi),  so  that  *  1  and  AT  «  A U^. 

Finally,  after  the  tensor  components  of  stress  were  found,  the 
physical  components  of  stress  were  calculated.  At  the  neck,  where 
shear  strain  does  not  occur,  the  dimensionless  physical  components  of 
stress  are  simply: 
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The  dimensionless  total  load 
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was  calculated  at  the  neck  by  intea 


grating  ( — — )  over  the  current  area  A  of  the  neck;  namely, 


2-  .D2  J  % 


(43) 
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Since  dA  =  2*( r  +  u  )d(r  ♦  u  )  =  2jt(1  +  f  +  r2h)(l  +  f  +  3r2h)rdr,  equation 

r  r 

(43)  becomes 
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with  f  and  h  being  evaluated  at  x  =  0.  The  value  of  the  total  load  was 
then  obtained  by  the  numerical  integration  of  (44). 

VI.  Results  and  Discussions 

The  results  presented  are  based  on  the  following  conditions.  The 

initial  half-length  of  the  specimen  is  L  =■  2,  and  the  initial  radius  is 
2 

R(s)  *  1  -  ce  *  t  with  the  imperfection  arbitrarily  chosen  as  c  -  0.005 

(the  results  are  not  sensitive  to  the  magnitude  of  «).  The  Strain¬ 
s'  , 

hardening  exponent  is  n  =  8,  the  ratio  (-^-)  i*  0.0072  (E  =  25.7  x  10  psi. 
and  9  =  .186  x  10^  psi. ),  and  Poisson's  ratio  is  v  =  -j.  The  results  in¬ 
clude  the  total  load,  the  displacements  and  the  stress  and  strain  distri¬ 
butions  calculated  at  every  stage  of  deformation.  Each  stage  of  deforma¬ 
tion  was  specified  by  a  different  value  of  the  engineering  axial  strain 
UL  UL 

-j— .  In  the  present  paper  -j—  was  varied  with  increment  .0025,  from 
0  to  .  50,  at  which  point  the  calculation  was  terminated  for  reasons 
which  will  be  given  later.  The  spatial  mesh  sizes  used  in  the  numerical 
analysis  were  =  <05  and  =  .025. 

6. 1.  The  total  load 

j 

Fig.  8  shows  a  plot  of  the  dimensionless  total  load  — j —  (as  cal- 

UL  ***0*8 

culated  by  (44))  against  the  engineering  strain  -s— .  The  total  load  reaches 
UL 

its  maximum  at  -7—  =  0.  145.  It  is  interesting  to  notice  that  the  plastic 

UL 

unloading  (which  starts  at  the  ends  of  the  bar)  does  not  begin  until  -j— 

reaches  0.  180,  which  is  slightly  after  the  total  load  has  started  decreasing. 

T  UL 

For  comparison,  the  plot  of  —  •j —  versus  —j—  as  shown  in  Fig.  4 

'V. 

for  the  perfect  bar  was  reproduced  in  Fig.  8.  It  is  seen  that  the  maximum 
value  of  the  total  load  of  the  imperfect  bar  agrees  quite  well  with  that 
predicted  for  the  perfect  bar.  The  discrepancy  between  the  two  total 
loads  begins  shortly  after  the  total  loads  start  decreasing. 


FIG-  8  TOTAL  LOAD  VS-  ENGINEERING  AXIAL  STRAIN 
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The  relation  between  -  — y  —  and  the  radius  at  the  neck  R  is  shown 

*Ro'. 

in  Fig.  9.  The  results  for  the  perfect  bar  were  also  reproduced  from 
Fig.  5  for  comparison. 

6.  2.  The  displacements  and  the  deformations 

Fig.  10  shows  the  radial  displacement  up  on  the  lateral  surface, 
r  -  R(z),  at  different  stages  of  deformation.  During  the  earlier  stages 
of  deformation,  the  inward  radial  displacement  increases  quite  uniformly 
throughout  the  length.  But  for  later  stages,  the  inward  radial  displace¬ 
ment  increases  only  locally  at  the  middle.  At  the  end  of  the  specimen 
the  radius  expands  slightly  when  the  total  load  decreases  (sec  the  enlarge¬ 
ment  in  Fig.  10). 

The  axial  displacement  u  is  shown  in  Fig.  11  for  different  stages 
of  deformation.  In  this  figure  the  solid  lines  indicate  the  axial  displace¬ 
ment  on  the  lateral  surface,  r  =  R(z),  and  the  dashed  lines  show’  the  axial 
displacement  along  the  center  line,  r  =  0.  The  slopes  of  these  curves 

uL 

become  much  steeper  at  the  middle  for  large  values  of  -j—  .  This 
implies  that  the  axial  strain,  which  is  approximately  equal  to  the  slope 
of  these  curves,  becomes  concentrated  locally  at  the  middle  when  the 
deformation  is  large. 

The  localization  of  the  displacements  at  the  middle  of  the  specimen 
as  shown  in  Fig.  10  and  Fig.  11  can  easily  be  explained  qualitatively  as 
follows.  The  radius  at  the  middle  is  initially  smaller  than  the  radius  at 
the  ends  due  to  the  initial  imperfection.  As  the  elongation  proceeds,  the 
stress  and  strain  at  the  middle  cross-section  are  larger  than  at  the  end 
cross-sections;  therefore,  the  discrepancy  between  the  radii  increases. 
When  the  rate  of  decrease  of  cross-sectional  area  exceeds  the  rate  of 


FIG.  9  TOTAL  LOAD  VS.  NECK  RADIUS 


FIG.  10  RADIAL  DISPLACEMENT  ON  LATERAL  SURFACE 
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increase  of  the  average  axial  stresses  at  the  middle,  the  total  load 
begins  to  decrease.  Shortly  after  this,  the  effective  stress  decreases 
at  the  ends,  and  hence  plastic  unloading  begins.  Once  the  plastic  un¬ 
loading  occurs,  the  area  contraction  at  the  ends  stops  (in  fact,  the  radii 
at  the  ends  expand  slightly  as  a  result  of  the  decrease  of  the  axial  stress); 
therefore,  the  discrepancy  between  the  radii  of  the  middle  and  the  ends 
accelerates.  As  the  elongation  increases  further,  the  stress  at  the  neck 
becomes  larger,  but  the  cross-section  contracts  at  an  even  faster  rate, 
so  that  the  total  load  decreases  further.  This  results  in  larger  plastic 
unloading  regions  at  the  ends  and  aggravates  the  discrepancy  in  the  radii; 
thus,  the  displacements  become  localized  at  the  middle. 

The  location  of  the  boundary  between  the  plastic  unloading  region 
and  the  plastic  loading  region  is  shown  in  Fig.  12  at  different  stages  of 
deformation.  To  the  left  of  this  boundary  is  the  plastic  loading  region  and 
to  the  right  of  this  boundary  is  the  plastic  unloading  region.  This  boun¬ 
dary  moves  toward  the  middle  as  the  total  elongation  increases. 

From  the  knowledge  of  the  displacements,  the  shapes  of  the  necked- 
down  specimen  at  different  stages  of  deformation  were  calculated  and  are 

U, 

shown  in  Fig.  13.  At  -r-  =  .  50,  the  localized  necking  results  in  a  53% 

Lj 

reduction  of  the  radius  at  the  middle,  while  the  reduction  of  the  radius 
at  the  end  is  only  about  6.  3%.  The  detailed  deformation  in  the  specimen 
can  be  seen  from  Figs.  14  (a),  (b)  and  (c).  These  figures  show  the  de- 

UL 

formations  of  a  square  grid  of  the  original  body  (see  Fig.  7)  at  — —  =  0.  30, 

Lj 

0.  40  and  0.  50. 

The  shapes  of  the  necked  profiles  of  the  tension  specimen  shown  ob¬ 
viously  resemble  the  shapes  observed  in  tests.  However,  no  precise 
comparison  has  been  made. 


=0.0072 


UNLOADING  REGION 


» 


-42- 


O 

z 

Q 

< 


Z 

o 

o 

Ul 

(t 


h- 

CO 

< 


o 

z 


p 


o 
a  ui 
Z  a: 


co 

< 

5? 


4*1 


i 


O 


TTWT* 


f 

! 


-45- 

6.3.  The  strain  and  stress  distributions 

The  stress  and  strain  distributions  have  been  calculated  in  the 
whole  specimen  during  its  loading  history.  Only  some  of  these  dis¬ 
tributions  at  the  neck  will  be  shown. 

Figs.  15,  16  and  17  show  the  Lagrangian  strain  tensor  components 

UL 

across  the  neck  section  at  -j—  =  .  25,  .  30,  and  .  35.  The  maximum  axial 

strain  (r?  -)  is  seen  to  occur  at  the  center.  It  is  interesting  to  recall 
j  -3  m  e  x 

that  the  Bridgman  [1]  and  Davidenkov  and  Spiridonova  [2]  analyses  were 
based  on  the  assumption  that  the  axial  strain  at  the  neck  was  uniform. 

In  the  present  analysis  after  the  Lagrangian  strain  components 

are  found,  the  value  of  *\  —  ,  which  was  assumed  to  be  close  to  1  in  order 

'  g 

i  j  ^  [  j 

to  use  the  approximation  <r  J  =  q  J,  was  calculated.  The  maximum  error 

UL 

of  this  assumption  was  found  to  be  less  than  1.  5%  for  — —  under  .  50. 

The  distributions  of  the  physical  components  of  stress,  °’rr.)  '’’qqj 

and  <r  ,  at  the  neck  are  shown  in  Figs.  18a,  19a,  and  20a  for 
UL  ZZ 

—  =  .  25,  .  30  and  .  35,  and  the  distribution  of  the  effective  stress  at 

these  stages  of  deformation  are  shown  in  Fig.  18b,  19b,  and  20b.  The 

UL 

analysis  showed  that  for  -7 —  less  than  .27,  (<r  )  occurs  at  the 

7  L  tt  zz  max 

center  as  shown  in  Fig.  18a,  and  for  -s —  over  .27,  (<r  )  moves 

between  the  center  and  the  peripheral  surface  as  shown  in  Figs.  19a  and 
20a.  The  stress  distributions  shown  may  be  inaccurate  at  comparatively 
large  deformations,  since  then  <r  does  not  vanish  on  the  peripheral  sur¬ 
face  at  the  neck.  (The  requirement  c t  =  0  on  the  peripheral  surface  at 

the  neck  is  a  natural  boundary  condition  of  the  variational  principle;  to 
impose  this  condition  is  beyond  the  control  of  the  present  calculation  pro¬ 
cess.  The  error  in  c  on  the  surface  can  be  used  as  a  criterion  for 

r  r 

checking  the  accuracy;  this  accuracy  can  be  improved  by  including  more 
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FIG.  15  LAGRANGIAN  STRAIN  COMPONENTS 
AT  THE  NECK  FOR  -^=.25 
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FIG.  16  LAGRANGIAN  STRAIN  COMPONENTS 
AT  THE  NECK  FOR  ^  8-30 
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FIG.  17  LAGRANGIAN  STRAIN  COMPONE 
AT  THE  NECK  FOR  ^=.35 


FIG.  18A  STRESS  DISTRIBUTIONS  AT  THE  NECK 
FOR  UL  =,25 

T~ 


« 


-50 


FIG.  18  B  EFFECTIVE  STRESS  DISTRIBUTION 
AT  THE  NECK  FOR  UL  =,25 
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FIG.  19o  STRESS  DISTRIBUTIONS  AT  THE  NECK 
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FIG.  190  EFFECTIVE  STRESS  DISTRIBUTION  AT 
THE  NECK  FOR  ^  *  .30 


FIG-  20A  STRESS  DISTRIBUTIONS  AT  THE  NECK 
FOR  » .33 
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terms  in  the  power  expansions  in  r  in  the  assumed  displacements  (34) 
and  (35).  )  The  results  of  the  present  analysis  shcxv  that  «r  becomes 
negative  on  the  peripheral  surface.  However,  these  negative  values 
for  cr  are  localized  near  the  peripheral  surface  only.  Therefore  the 
stress  distributions  obtained  may  not  be  too  incorrect  inside  the  necked 


uL 

section.  At  — —  =  .  50,  the  negative  value  of  cr  on  the  surface  becomes 
L  ’  °  rr 

about  half  of  (cr  )  ,  and  for  this  reason  the  calculation  was  terminated, 

zz  max’ 

Also  shown  in  Figs.  18,  19  and  20  are  the  stress  distributions 
calculated  by  the  formulas  of  Bridgman,  and  Davidenkov  and  Spiridonova, 
on  the  basis  of  the  neck  shapes  predicted  by  the  present  analysis.  These 
stress  distributions  are  given  by  the  formulas 
(1)  Bridgman 

O’  rr- 


7Tcr  (R2  +  2Rp)  ln(l  +  4“-) 
s  £  p 


°rr  °ee  =  (  %  )  ln  R2  +  2Rp  -  r2 
°"s  ‘’’s  %  2Rp 


tr  cr  cr 
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(2)  Davidenkov  and  Spiridonova 
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Here  the  total  load  T  is  calculated  by  (44),  r  is  the  distance  frorr.  the 
center,  R  is  the  current  neck  radius,  and  p  is  the  radius  of  curvature 
at  the  neck.  Both  R  and  p  were  measured  from  the  calculated  shapes 
of  the  tension  specimen  shown  in  Fig.  13.  The  effective  stress  cal¬ 
culated  by  [1]  and  [2]  are  seen  to  be  close  to  the  effective  stress  calcu¬ 
lated  by  the  present  analys;s.  But  this  seems  to  be  a  coincidence  since 
the  detailed  stress  and  strain  distributions  of  the  present  analysis  were 
not  close  to  those  of  [1]  and  [2].  Perhaps  it  should  be  emphasized  here 
that  the  above  comparison  may  not  be  too  meaningful,  since  the  radius 
R  and  the  radius  of  curvature  p  at  the  neck  which  were  required  for  the 
calculation  of  the  stress  distributions  given  by  Refs.  [1]  and  [2],  were 
provided  by  the  present  analysis. 

VII.  Concluding  Remarks 

The  most  significant  result  of  the  present  analysis  is  that,  for  the 
first  time,  the  shapes  and  the  deformations  of  a  necked-down  tension 
specimen  were  rationally  calculated.  The  stress  and  strain  distributions 
in  the  specimen  duting  its  loading  history  were  also  found,  although  the 
numerical  results  for  them  are  still  far  from  conclusive. 

Some  improvements  of  the  numerical  technique  will  be  needed  for 
the  present  analysis  to  obtain  more  accurate  results.  Then,  on  the 
basis  of  the  present  analysis,  the  stress  and  the  deformation  of  a  tension 
specimen  throughout  its  loading  history  can  be  predicted  if  the  elastic 
constants  and  the  strain-hardening  characteristics  of  a  material  are 
known.  This  will  provide  knowledge  of  the  stress  conditions  in  a  necked- 
down  tension  specimen  before  it  fractures;  thus  it  perhaps  will  lead  to  a 
better  understanding  of  the  fracture  phenomena  in  a  tension  bar,  such 
as  the  cup-cone  type  fracture. 


APPENDIX 

Inversion  of  the  Constitutive  Equation 


Since 
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equation  (10)a  can  be  rewritten  as 
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We  multiply  both  sides  of  (45)  by  G  G  J  and  rearrange  it  to  get 
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On  the  right  hand  side  &e  and  cr^  are  unknowns.  In  order  to  find  ir^, 

we  multiply  (45)  by  G1*^  to  get 

Gijr?ij  =  ^  [(1  +  v)(crklGj<1  +  2<rklrjkl)  -  3v«rk] 
in  which  (crklGkl  +  2<rklnkl)  =  =  '  Rearranging  (47)  gives 


f 


-58- 


.  k 
rk* 


,■  “y  GlJh.. 

1  -  2v  ''ij 


(48) 


Next,  in  order  to  find  we  multiply  (10)a  by  s1^  and  obtain 
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With  the  use  of  the  relations  s!  =  <r!  -  4-6|<rk  s!  =  fr!  -  and 
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Note  that  in  (50),  (—  -  )  is  always  positive,  so  that  the  loading  condi¬ 
tion  ore  ^  0  can  be  replaced  by  ^  0.  Substitution  of  (48)  and  (50) 

into  (46)  gives  the  final  form  of  the  inverted  constitutive  relations 
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The  necking  process  of  an  axisymmetric  tension  specimen jg a^e  of  elastic  - 
plastic,  strain-hardening  material  is  analyzed  by  a  generalized  Jt- flow  theory  for 
large  deformations.  The  governing  equations  are  solved  in  a  Kamtbrovich  type 
approach  based  on  a  variational  principle.  The  effects  of  both  geometric  nonlineari¬ 
ties  resulting  from  large  deformations  and  the  physical  nonlinearities  arising  from 
plastic  material  behavior  are  considered.  Numerical  results  have  been  found  for  the 
stress  and  deformation  histories  in  the  specimen  up  to  a  50%  reduction  of  the  neck 
radius.  The  shape  of  the  neck,  for  the  first  time,  is  rationally  calculated. 
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